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Abstract 

This paper investigates the possible cosmological imphcations of the pres- 
ence of an antisymmetric tensor field 9 related to a lack of commutativity of 
spacetime coordinates at the Planck era. For this purpose, is promoted to 
a dynamical variable, inspired by tensor formalism. By working to quadratic 
order in 9, we study the field equations in a Bianchi I universe in two models: 
an antisymmetric tensor plus scalar field coupled to gravity, or a cosmologi- 
cal constant and a free massless antisymmetric tensor. In the first scenario, 
numerical integration shows that, in the very early universe, the effects of the 
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antisymmetric tensor can prevail on the scalar field, while at late times the 
former approaches zero and the latter drives the isotropization of the universe. 
In the second model, an approximate solution is obtained of a nonlinear ordi- 
nary differential equation which shows how the mean Hubble parameter and 
the difference between longitudinal and orthogonal Hubble parameter evolve 
in the early universe. 
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Typeset using REVT^ 



I. INTRODUCTION 



Spacetime noncommutativity is one of the key new hints which follow from recent de- 
velopments in quantum field theory. It has been recently realized, in particular, that a 
consequence of string theory [1,2] is that the structure of spacetime becomes noncommuta- 
tive [3], which can be described loosely as an analog of a quantum phase space, in terms 
of the algebra generated by noncommuting coordinates [a;'^,a;'^] = i9^^'^ with 9^^'^ an anti- 
symmetric tensor^. The idea behind spacetime noncommutativity is very much inspired by 
quantum mechanics. A quantum spacetime is defined by replacing canonical variables with 
self-adjoint operators which obey Heisenberg-like commutation relations [x^, x'^] = iO'^'^, and 
can be viewed as the smearing out of a classical manifold, with the notion of a point re- 
placed with that of a Planck cell. It was von Neumann who first attempted to rigorously 
describe such a quantum "space" and he called this study "pointless geometry", referring 
to the fact that the notion of a point in a quantum phase space is meaningless because of 
the Heisenberg uncertainty principle of quantum mechanics. This led to the theory of von 
Neumann algebras and was essentially the birth of noncommutative geometry, referring to 
the study of topological spaces whose commutative C*-algebras of functions are replaced by 
noncommutative algebras [3-8]. The idea of noncommutative geometry was revived in the 
eighties by Connes [5] and others, who generahzed the notion of a differential structure to 
the noncommutative setting, i.e. to arbitrary C*-algebras. A theory on a noncommutative 
space replaces the noncommutativity of operators associated to spacetime coordinates with 
a deformation of the algebra of functions defined on spacetime. In this context classical 
general relativity would break down at the Planck scale because spacetime would no longer 
be described by a differentiable manifold, and at these length scales quantum gravitational 
fluctuations become large and cannot be ignored. We stress, however, that the form of non- 
commutative geometry we are interested in is not directly related to current string theories 



Hereafter we use the natural units % = c = 1. 
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(see below). 

In the past few years several authors [9-11], including some of the present authors [12-16], 
have considered the possible effects of noncommutative geometry and Planck scale physics 
in cosmology. In particular, it has been shown that deformation of spacetime and/or phase 
space algebras may lead to several interesting features in the power spectrum of primordial 
perturbations produced during the inflationary era [17-20]. In all these investigations, how- 
ever, 9^^ has been taken to be constant, or with an a-priori modelled time evolution. In 
view of general covariance one may expect that d^'' should be rather considered as a dy- 
namical tensor, coupled to gravity and possibly affecting the cosmological evolution of the 
early universe. This is actually a crucial point which deserves a thorough treatment, here 
summarized by relying in part upon Ref. [16]. On the one hand, it is true that, if one looks 
at the interplay between string theory and noncommutative geometry, one has to consider a 
constant 5-field and hence a constant ^^'^ [21]. On the other hand, the commutator of the 
is a tensor, whose transformation under boosts yields nonvanishing spacetime components. 
Moreover, possible violations of unitarity are taken care of by imposing the conditions [16] 

Prom the point of view of general formalism, it is therefore legitimate to address the question 
of whether a broader noncommutative picture can be consistently built. This implies a 
departure from current models which only exploit a constant 9^" , and suggests starting 
from a nonlocal action functional with * product of fields in the presence of nonvanishing 
spacetime curvature [16]. We therefore assume, hereafter, that the background geometry 
remains a classical pseudo-Riemannian geometry endowed with a Levi-Civita connection V, 
while the * product of scalar fields at the same spacetime point is defined by [16] 

OO -1 

^{x)*^{x) ^ ^-(i/2)'=r^'^H^)-^'^'='^H^)(V;.,...V^,^)(V.,...V.,V') ■ (1-1) 

fc=0 

Similarly, having treated classically the geometry, we assume for tensor fields that 

Fm..Xs * F^'-^' = /^'^^../^'■^Fa,..A. * Fr,...r^ , (1.2) 
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where 

oo 1 
fe=0 

As is stressed in Ref. [16], the occurrence of covariant derivatives in our definitions (1.1) and 
(1.3) spoils associativity of the * product. However, noncommutative effects are already 
present at quadratic order in 9^^ , and our O^'' will be taken to be sufficiently small so that 
higher order terms in the action functional are negligible. 

Another important goal to be pursued is a fully consistent study of primordial per- 
turbations, which should take into account the 6^^ perturbations as well, generalizing the 
formalism of gauge-invariant perturbations [22] to anisotropic background metrics. 

As a first step in this programme, in this paper we discuss the possible dynamical evolu- 
tion of a background, time-dependent antisymmetric tensor [23-30] [31-38] in two possible 
scenarios: a free massless ^^'^ in presence of a cosmological constant, the latter being intro- 
duced as the easiest way to trigger an infiationary dynamics, and a more general scenario 
where the noncommutative field 9^^ is coupled to a scalar inflaton. As mentioned the 
presence of 9^^ breaks the isotropy of the universe, and hence only spatial homogeneity 
is preserved, leading in turn to a dependence of 9^" on time only. In this framework the 
appropriate geometry for the universe is therefore a Bianchi I model. In Sec. II we obtain 
a nonlinear system of the background equations in the presence of infiaton plus a coupling 
term between infiaton and 9^" , and the energy-momentum tensor. In Sec. Ill we consider 
the simpler cosmological term model which however can be worked out analytically, at least 
for those initial conditions which are of some interest and may lead to an early stage where 
the energy-momentum tensor is dominated by 9^^ . Concluding remarks and open problems 
are presented in Sec. IV, while relevant details are given in the Appendix. 
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II. ANTISYMMETRIC TENSOR PLUS SCALAR FIELD COUPLED TO 

GRAVITY 



In order to describe a field dynamics which might lead to anisotropy in the early universe, 
the appropriate model is a Bianchi I universe (as we stated before) if one wants to preserve 
the spatial homogeneity. In this case the line element can be written as 

ds'' -Y,a^,{t){dx'f , (2.1) 

i=l 

and correspondingly the nonvanishing connection coefficents are (no summation over i is 
here meant) 

r°. = 5i>,a, , roi = - , Vz,j = 1,2,3 , (2.2) 
and for the Ricci tensor one has 

K=-j:'f , Bi = -si(^+^j:^ ■ (2-3) 

i=l Y'i k^i «fe / 

We consider a model in which there are both the antisymmetric tensor responsible of non- 
commutativity of spacetime and a minimally coupled massive scalar field which drives the 
inflation. The corresponding nonlocal action reads 



j d^xy/^ 



2 



R 1 1 .,, m 



16nG 12 ■ 2'^''' 2 

(2.4) 



where the part involving A is here introduced to mimic a 'time-dependent' mass term for 
O"", and 

H^^^a = + ^vOa^ + V^^^. (2.5) 

is the field strength associated to the antisymmetric tensor Of^i, — —9un (hereafter, Greek 
indices run from through 3, whereas Latin indices run from 1 through 3). It should be 
noticed that the kinetic term for O^^i, is inspired by a generalization of the Maxwell theory 
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[34] and that in Eq.(2.5) only the effects of partial derivatives survive ^. Now, by virtue of 
the definitions (1.1)-(1.3), one finds 



(2.6) 



(2.7) 



(2.8) 



Thus, to second order in 9^'^, which is relevant since 9'^^ is taken to be sufficiently small, 
only kinetic and potential term for the scalar field contribute, but with vanishing coefficient, 
since the former changes by the amount [16] (integration by parts yields also a third term 
which however vanishes if <^ = ip{t) only) 



and the latter changes by the amount [16] 

2 

^^rn^^J d'x^g 9^-'9^-Rl^^{dr^){d,^) . 



(2.9) 



(2.10) 



Since depends only on the time variable, both (2.9) and (2.10) vanish in our Bianchi I 
background. Thus, to quadratic order in 9^"^, we end up with the local action functional 

R 1 „ 1 A 



S = j d'^x^/^ 



(2.11) 



where V{(p) = hereafter. At this stage, the resulting energy-momentum tensor is given 



by 



(2.12) 



^Our notation agrees with the one used, for example, in Ref. [35]. 
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By using the expressions (2.11), (2.12) the resulting equations of motion are 

s^''R^,^ + 2X^%^ = , (2.13) 

VmV'^¥' + A<^V^'^'^ + ^ = , (2.14) 

R^,u-lg^..R = SttCT;, . (2.15) 

To be consistent with the spacetime homogeneity ansatz we assume that all fields are de- 
pending on time only. In this case the above equations read (hereafter Hi = hi/ai) 

A ^0^ = 0, (2.16) 

'^i3 + (E ^fc) % - 2 (i/i + Hj) % + 2 A % = , (2.17) 

V3+ [j^H^^ + Xipe^^e^^ + ^-^^Q , (2.18) 

3 

E {^^k Ojk - 2 A g,^ ^^.^^ _ Q ^ ^-^j ^2.19) 



k=l 



E ^ = 8 TT G (vi^) -^'-^V' O.J'"'^ , (2.20) 

- + -E- = 8^G 1 jj jj^.. + y _ ^ ^^M^ _ 1 H,.^ff- + 2 A e.^r'l . (2.21) 
di a-i t^^ak \6 2 2 y 

Of course, Eqs. (2.16) and (2.17) result from (2.13), whereas (2.18) is obtained from (2.14). 

The remaining equations are the Einstein equations where, in particular, (2.19) provides a 

Bianchi I universe. 

Using Eq. (2.16) one easily gets 9^^ — 9oi — 0. By virtue of Eq. (2.19) one can show that 
the only possible solution has only one nonvanishing component of 9ij, e.g. 9i2- Moreover, 
since isotropy is broken and the residual invariance is 50(2), it is rather natural to choose 
ai — a2 = a±, = ttL, with corresponding Hubble parameters H± = a±/a±, Hl = ciL/aL- 
In this case the equations become 

9i2 + {Hl - 2Hj_) + 2 A ^12 = , (2.22) 
^ + {Hl + 2H^) + A^ + ^ = , (2.23) 



'1 .2 , , ^9^ 



HI + 2H^Hl^^'kG \ :-^' + V{^) + :-"-^ + X^^\ , (2.24) 



H^ + H^{HL + 2H^_) = 87rG {V{cp)+ X . 
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Hr+Hr(Hr+2H 



8nG (^'^ 
\ a 



/o2 Q'2 
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(2.25) 
(2.26) 



Let us define the mass parameter /j, = mpi/\/%T^ — l/VSnG. In terms of this quantity we 
can write a dimensionless system of differential equations. By defining 

X = (p/fji, y = 9i2/{ijLa\), H± = H±/ijl, Hl = Hl/h, V{x) = V {/ix) / /i'^ , 



and using the dimensionless time t = t/i we get 

/ + (Hl + 2i?x) y+2 (Ax^ -2Hl + V{x) + Xx'^y^) y = 

AT/ 

" + (Hl + 2H^) x' + Xxy^ + — ^0 , 

2~ ' ■ ' 2 
H'^ + H± (Hl + 2H^ = V{x) + Ax' 



X 



/2 1 

Hi + 2HJIl = ^ + V{x) + ^ (y' + 2H^y)' + A 



2 2 

X y 



H'l + Hl {Hl + 2H^) = (y' + 2i/^t/)' + V{x) - X 



2 2 
X 7/ 



(2.27) 
(2.28) 

(2.29) 
(2.30) 
(2.31) 



where the 'prime' denotes the derivative with respect to r. The truly independent equations 
are given by (2.27), (2.28), (2.30) jointly with 

- I +V{x) + ^ {y' + 2H^ 

2H^ 

On recalling the e-folding definition in the orthogonal direction, i.e. N±{t) 
log(aj_(T)/a_L(ri)), one can easily prove that 



^Hl + 2H^)=^(^^ +V{x) + ^ {y' + 2H^yy +Xx'y' + 3Hj 



(2.32) 



— ^H 

dT ^ dN^ ' 

^ tt2 d'' ~, d 
dr'^ ^ dNl ^ ^ dN^ ■ 

By using N±_ as the evolution parameter and defining z = H^ we find 
d-'y 



dN^ 
d^ x 



2 2\ dy 



^ dNl 
dz 



(^ri \ ^ 2 2\ dx ^ 1 
+ \y{x) + Ax'y'j 3^ + Ax?/' + 



dN^ 



= T/(x)+Ax2y' 



dN\ 



+ 2(Ax'-2z + T/(x) + A 
, 

+ 2?/) +6 



a;2 y^\y = ^ 



dx 



dN 



+ 



5x 

' dy 
dNl 



(2.33) 
(2.34) 

(2.35) 
(2.36) 
(2.37) 



In the following figures we show, on choosing different initial conditions, that a range of A^j_ 
exists in which the 9i2 component of the 9 field dominates on the ip field; of course, this 
results from a special choice of initial conditions (see following section for a more thorough 
discussion of underlying issues). Moreover, Hl dominates as well, i.e., there is anisotropy. 
Outside this range, 612 approaches zero and (p dominates driving the isotropization, which 
can be seen in the figures where Hl and H± reach the same constant value. 
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FIGURES 
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FIG. 1. The X-axis corresponds to the e- folding parameter N± and on the y-axis we plot if^ 
and H± with dashed and continuous line, respectively; we take V{x) = qx'^, with initial conditions: 
Hi = 10, 9i2 = 100, ^ = 0, (/; = 10, ^ = 0, and coupling constants A = 10"^, q = 10. 
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o l , -^rrr^r-^, 

5 10 15 20 25 

FIG. 2. The x-axis corresponds to the e-folding parameter N±, and on the y-axis we plot 612 
and (f with dashed and continuous line, respectively, with same initial conditions of the previous 
figure. 



FIG. 3. We plot Hl and H± as functions of the e-folding parameter N±, with V{x) = qx^, 
and initial conditions: Hj_ = 10, ^12 = 100,^^ = l,ip = 10,^^ = 0, and coupling constants 



A = -10-5,g = 1. 
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FIG. 4. The fields 612 and tp are plotted as functions of the e-folding parameter N±. 



III. COSMOLOGICAL CONSTANT AND MASSLESS ANTISYMMETRIC 

TENSOR 

While the equations in the previous section can only be solved by means of numerical 
methods, one can envisage a simpler scenario which embodies the main features but allows 
for an analytic approach. For this purpose we consider a Lagrangian density for the anti- 
symmetric tensor and gravity including a cosmological constant term (with Hf^i^p defined 
in Eq. (2.5)) 



IGnG 
1 



IQttG 



12 



+ 0(^3 



(3.1) 



Thus, on working to quadratic order in 6"*'^, the resulting action is invariant under the gauge 
transformations 



(3.2) 



and the 9^^ dynamics can be also recast in a different form by introducing the pseudo-scalar 
Kalb-Ramond field x via the duality transformation d^x = '^iiupaH^'^" , with e^yfj„ the volume 
four-form. Actually, since we are assuming homogeneity and a purely time-dependent 9 field 
this implies a time-independent x fi^ld with a linear dependence on spatial coordinates (see 
Eq. (3.15)). 
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Using the results of the previous section in the particular case we are considering, the 
resulting equations of motion are 

v^a^ua = , (3.3) 
R^^u-lg^^uR = 8TiGT^,{9) + Ag^, , (3.4) 

and the energy-momentum tensor reads 

THie) = \h,^,H''^^ - ^^H.^.H"^" . (3.5) 

As in the previous section, we assume Bianchi I geometry. The equations of motion (3.3) 
and (3.4) now read 



% + [12 Hkj Oij -2{Hi + Hj) % = , (3.6) 

i2^ikejk = , Vt^j (3.7) 
fe=i 

E^ = A, (3.8) 

i=l 

h + ^J2^ = A-4nG (h,..H^- - I H^.^H'^-) . (3.9) 

Unlike the previous section, Eq. (2.16) is now replaced by a gauge- fixing condition. In the 
following we will use the gauge 9oi = 0. A detailed analysis of the admissibility of this gauge 
fixing is reported in the Appendix. 

The effect of (3.6) and (3.7) is to restrict severely the possible choice of initial conditions 
for 9ij. In particular either all ^ij(O) values vanish and thus 9ij is a constant tensor, or at 
most only one component %(0) is nonvanishing. Since the antisymmetric tensor enters the 
field equations only through its first- and second-order time derivatives the case of a constant 
^-tensor reduces to a pure A-term cosmology and is thus uninteresting for our analysis. We 
therefore consider the case ^13(0) = ^23(0) = 0, ^12(0) 7^ 0, and we eventually get 

O12 + {Hl - 2H^_) 012^0 , (3.10) 



2 

Hl + 2H^HL = A + 47rG"-f , (3.11) 

a I 
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Hi_ + {Hl + 2H^) = A , (3.12) 



o2 

HL + HL{HL + 2H^)=A + 87rG^ . (3.13) 

a I 



Equation (3.10) determines the behavior of 6*12 in terms of the scale factors 

012 = ^12(0)^ , (3.14) 

where in (3.14), by virtue of the arbitrariness in the choice of 0^(0) and a±{0), we have 
assumed, without loss of generality, 0^(0) = aL(0) = 1. Notice that this immediately gives 
the spacetime dependence of the Kalb-Ramond scalar field 

X^0i2(0)x' + xo , (3.15) 

where, according to (2.1), denotes the longitudinal spatial coordinate. The equation 
(3.15) is actually consistent with spatial homogeneity, while it breaks isotropy. 

The previous equations can be recast in a simpler and more useful form by introducing 
the variable 

e = ^log(aiai) , (3.16) 
which represents the average e-folding. Hence one gets 



c 



.2 



HI + 2H^Hl^A+^ , (3.17) 

"l 



2a 



{H'^ + 3H^) {Hl + 2H^) = 3A , (3.18) 
{H'l + 3Hl) {Hl + 2Hj_) = 3A + 3— , (3.19) 



where = 8nG (^12(0)) and ' = d/d^. It is also convenient to set 



and hence from (3.16) 



MO=exp|e+^| , (3.20) 



a±(0=exp|e-^} . (3.21) 
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Furthermore, we introduce the mean Hubble parameter, H = [Hl + 2H±)/3, and the 
asymmetry function h = — H±. Since 0^(0) = a±{0) = 1, Eqs. (3.20) and (3.21) yield 
fl{0) = 0. Thus one gets 

^ (h')' + 3H' = A + j exp {-2^ - m} , (3.22) 
/i' + 3/i = ^exp{-2e-Q(0} , (3.23) 
H'-^h'^^ + ^e^p{-2^-n{0} , (3.24) 

where the last equation is the Hamiltonian constraint. Actually the functional dependence 
of h in terms of H and Q can be already obtained by its very definition 

h = H(^-^]^ -Hn' , (3.25) 

VOL a±/ 4 

showing that Q! is directly related to the isotropy breaking. In the following we will consider 
an expanding universe in both longitudinal and transverse directions, Hl, H± > 0. In this 
case we note that —2 < Q' < 4. We will use this information later on. Notice that fl' = 4 
would correspond to H± — 0, which is forbidden by (3.17). 

The equation (3.22) is a first-order inhomogeneous equation whose solution reads 

= i^H^O) - ^) exp (-60 +lc' exp (-60 /^^ exp (4C - m) dC , (3.26) 

which, together with (3.25), gives the Hubble parameters in terms of the function 0(0- The 
latter can be determined as the solution of a second-order differential equation, obtained by 
inserting (3.25) into the differential equation (3.23) and the Hamiltonian constraint (3.24). 
The latter operation yields 

If the differential equation (3.22) is also exploited to express HH', we eventually find 



\ _ ^li^] exp (-2^-0(0) ^ 
, 16 y 2A + c2 exp (-2^-0(0) 
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with initial conditions Q(0) = 0, Q'(0) = (4/3)/i(0)///(0). 

This equation can be hardly solved analytically in the general case, but it reduces to 
much simpler forms in the two regimes, when either the cosmological constant or the O^i, 
fields dominate the energy-momentum tensor. At very early times the latter is likely to be 
largely the dominant component. We can rewrite the parameter (? as — mj^Q/m^pi, with 
TUNC the scale where the classical picture of spacetime manifold breaks down. The order 
of magnitude of A can be instead constrained by the fact that it drives infiation in the late 
stages, and represents in the slow-roll approximation the potential of the infiaton field. This 
is severely bounded by the fact that it should account for the correct amplitude of primordial 
perturbations. For example, for a polynomial potential = A0"/n, the requirement of 
slow-roll dynamics and perturbation amphtudes of the order of 10~^ gives A < 10~^^mp; 
[38]. As long as m]\fc ^ mpi the value of A is several orders of magnitude smaller than c^, 
so in this case the early dynamics is fully determined by Of^i,. Smaller values of ttinc, such 
that ~ A, cannot be ruled out of course, but in this case Eq. (3.27) can only be solved 
numerically. 

Hereafter we specialize to the case m^vc ~ it^pi- At early times therefore Eq. (3.27) 
takes the simplified form 

Q"(0 + 2(Q'(0 - 4) (l - ^) = . (3.28) 

Since the 9^1, contribution is diluted with expansion as a|, this equation holds approxima- 
tively for values of ^ smaller than the value ^* such that A ~ c^/a|(^*), i.e. 

e* + ^-log^. (3.29) 

Later expansion is instead driven by A and hence, by neglecting the 9^,, contribution, we 
have 

O"(e) + 3O'(0(l-^^) =0 . (3.30) 

This equation of course should give back the isotropization phase leading to a de Sitter 
phase. 
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We begin by studying the early time evolution. Before doing this it is worth discussing 
the values of initial conditions for Hubble parameters. Prom Eqs. (3.12) and (3.13) we 
see that 9f^i, acts as source for gl only, the evolution of a± being expected to be much 
slower. In other words, the antisymmetric tensor drives the expansion of the longitudinal 
scale factor only. The most natural choice at ^ is therefore Hl{0) >> H±{0), that is to 
say Q'{0) = 4 - 12e, with e = H±{0)/Hl{0) « 1. In this case the solution of Eq. (3.28) is 
particularly simple. On defining A(^) = Q(^) — 4^, the latter reduces to 

■ (3.31) 

Upon considering logA'(^) = y{^), this equation is solved exactly by separation of variables 
and subsequent integration, and yields (C being an integration constant) 

^ + C = - exp(-y) - I + 1 log (l + i exp(y)) . (3.32) 

On choosing the initial conditions 

A(0) = 0, A'(0) = -12e , (3.33) 

the approximate solution at small e reads, with a very good accuracy, 

A(e) = -^-^log(l + 12 6(l-^6)e) . (3.34) 

For example, we have checked numerically that this solution is accurate at the per thousand 
level up to ^ = 10^, if e < 10~^. These ranges fully cover the early stage. In fact Eq. (3.28) 
no longer holds at see (3.29), which, in view of the logarithmic behavior of A(^), for 
e < 10~^ is approximatively fixed by 

e*~^log^. (3.35) 

If c^/A < 10^^ we get < 15. On using (3.34) it is now possible to determine the evolution 
of the Hubble parameters as functions of ^. From Eqs. (3.26) and (3.25) we get 



H'iO ~ exp(-60 



(3.36) 



h{0 ~ 3exp(-3e)i/(0)[l - 3e] . (3.37) 
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Late-time evolution is ruled by Eq. (3.30). Since J1'(0 < 4 it therefore follows that 
asymptotically fl reaches a positive constant value D,^, so that h vanishes in the large-^ 
limit as expected. In particular, for large ^ Eq. (3.30) can be linearized, i.e. 

Q"(0 ~ -3n'(0 , (3.38) 
and hence it is immediate to get the approximate solution 

H'iO - f + {h'{Q - I) exp(-6(e - , (3.39) 
HO ~ KC*) exp(-3(e - e*)) . (3.40) 

It should be stressed that the action functional studied in the present Section or in Sec. 
II is not a low-energy limit of string theory with a constant dilaton, since the starting point 
remains a nonlocal action functional with * product of fields. Readers interested in string 
cosmology can be referred, for example, to the work in Refs. [39,40]. In particular, the 
work in Ref. [41] has studied string cosmology with a time-dependent antisymmetric tensor 
in a Bianchi 1 universe, but as is stressed in Sec. IV of Ref. [41], since a pure radiation 
plus dilaton solution has (f) — > const, a late-time isotropic radiation-dominated solution is 
a contracting universe. By contrast, in our model, when the universe reaches an isotropic 
stage it is still expanding. 

IV. CONCLUDING REMARKS 

Motivated by cosmology and noncommutative geometry, we have investigated the effects 
of an antisymmetric tensor in a Bianchi I early universe. What we do only holds at an 
intermediate stage where departure from ordinary spacetime geometry can be appreciated, 
while lack of associativity of the resulting * product in curved spacetime is negligible [16]. 
Moreover, such a stage implies a departure from the models with constant 9'^" which are 
relevant for string theory [21]. With this understanding, our original results are as follows. 

(i) In the first model, where the antisymmetric tensor 6 resulting from noncommutative 
geometry and a minimally coupled scalar field driving inflation both occur, we have found 
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by numerical methods that a suitable range of the e-folding exists such that prevails on the 
scalar field and the longitudinal Hubble parameter prevails on H^. This is the anisotropic 
era, but for larger values of the e-folding the antisymmetric tensor is damped and the scalar 
field dominates, leading in turn to isotropization. The very existence of the anisotropic era 
depends of course on the initial conditions chosen, but it appears interesting to have shown 
explicitly how the corresponding model can be built. 

(ii) In the second model, again in a Bianchi I universe, gravity with a cosmological constant 
A is coupled to an anti-symmetric tensor. The resulting nonlinear system of equations for 
the Hubble parameters H and h has been solved by first obtaining the differential equation 
(3.27) for the unknown function Q, and then finding approximate solutions when the effects 
of 6fj,i, prevail upon A, or the other way around. An accurate analytic description of the 
process leading to an isotropic final state of the universe has been therefore obtained. 

An outstanding problem in sight is now the development of cosmological perturbations' 
theory [38] within such a framework, with the hope of obtaining quantitative information 
on the effect of 9 on the formation of structure in the early universe. 
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APPENDIX 

In this Appendix we consider the admissibility of the gauge condition O^i — from the 
point of view of constraint analysis, and we begin with the simplest model, with action 
functional (as in Sees. II and III, the starting point is a nonlocal action functional, which 
reduces to a local action by retaining only quadratic terms in 9^") 
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J C^f^d^x , (Al) 

with Lagrangian density C = ^H^'^^'H^^p in a Bianchi I background, and field strength 
given by cycHc permutations of terms, according to Eq. (2.5). Only the effects of partial 
derivatives survive in and hence we find, by virtue of (2.1), 

C = ^H^'^'H,,, = -/'9'%k,^e^i,o , (A2) 

since the assumption of spatial homogeneity implies that 9'^''' can only depend on t. The 
term g^-'doifi^jfi is weighted with coefficient 3 — 3 = and hence does not occur in Eq. 
(A2). The Lagrangian L is obtained from C by means of (recall that the lapse function N 
is equal to 1 in a Bianchi I background, and that h denotes the determinant of the induced 
three-metric) 



J CVhd^x , (A3) 



and hence we find 



^'^^^ = ^9-9''0rs,o , (A4) 

7r°^^-^«.0. (A5) 

Equation (A5) deserves some comments: since the Lagrangian is independent of ^oi,o, the 
momentum conjugate to 6'oi vanishes. More precisely, the 7r°* can be seen as 3 primary 
constraints arising from the structure of the Lagrangian; as such, their vanishing is only a 
weak equation (?» 0), because 7r°' are well defined over the whole phase space of the theory, 
and only vanish on the constraint sub- manifold [42] . 
Our Lagrangian (A3) reads 

L = \ ^7r*^ %,oc?'x , (A6) 
with corresponding canonical Hamiltonian 
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I n'^ d'x-L = tt'^ %,o d'x . (A7) 



By virtue of the primary constraints (A5), however, the Hessian matrix is singular, and the 
time evolution is only well defined when the effective Hamiltonian H is considered. The 
latter is given by He plus a linear combination of primary constraints, i.e. [42] 



H 



I [l^''0^J,O + ^l%^)d'x . (A8) 



Moreover, we are still free to impose supplementary (more frequently called 'gauge') condi- 
tions, here chosen in the form 

Ooi^O . (A9) 

By doing so, we choose to regard the gauge conditions as constraint equations, in much the 
same way as the Coulomb gauge can be treated as a constraint equation in Maxwell theory 
[43] . We are therefore working with an extended Hamiltonian 



= J (^7r^^%,o + /^Voi + X%i^ d^x , (AlO) 



where the Lagrange multipliers ii\ A' can be evaluated by requiring preservation in time of 
the primary constraints ttoj and gauge constraints ^oi- For this purpose, note first that Eq. 
(A4) yields 

2 

%,o = -^T^ij , (All) 

and hence (with p = 

He = J (^7r^%,,- + /xVo, + A*^o^) d^x . (A12) 

All our constraints are then trivially preserved, without giving rise to further constraints, 
because 

^■Koi = {noi{x, t), He} ^ j X^{y, t) {noi{x, t), 9oj{y, t)} d^y = -Xi{x, t) , (A13) 



21 



d r 

-^o^ = {0o^{x, t), He] ^ J ^\y, t) {9o^{x, t), noi{y, t)}d'y = t) . (A14) 

Note that, by virtue of our gauge constraints (A9), the set of constraints has been turned 
into the second-class, a feature shared by all field theories after a gauge condition has been 
imposed [43], [44]. 

For the model studied in Sec. Ill, the full Hamiltonian constraint (3.11), when expressed 
in integral form, is the sum of (A12) and of the gravitational contribution. The latter is 
obtained from spatial integration of 

(167rG)G,,.,7r tt + + > 

where Gijki = ^^{Kkhji + huhjk — hijhki) is the DeWitt supermetric on the space of Rie- 
mannian geometries on S [45], tt'^ is the momentum conjugate to the induced three- metric 
and ^^^R is the three-dimensional scalar curvature (our sign for such a curvature is opposite 
to the one of Ref. [45]). 
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